<html>
<head>
     <title>Help: wsedecomp</title>
</head>
<body>
     <center><table><tbody>
          <tr>
             <td style="border: 1px solid rgb(170, 170, 170); padding: 5px; background-color: #ddd;">
                  <h4>WSEdecomp</h4>Split a signal into its downward/upward propagating and stationnary components
                  <br><small>Last modified: 25-Sep-2009 16:28:42</small>

          </tr>
          <tr>
             <td style="border: 0px solid rgb(170, 170, 170); padding: 5px; background-color: #ddd;">
                  <a href="http://guillaumemaze.googlecode.com/svn/trunk/matlab/codes/stats/WSEdecomp.m">Download here</a>
             </td>
          </tr>
          <tr>
             <td style="border: 1px solid rgb(170, 170, 170); padding: 5px; background-color: #fff; font-size:12px;">
                  <br>% WSEDECOMP Split a signal into its downward/upward propagating and stationnary components<br>%<br>% [Sf A FX FT] = WSEdecomp(S,DT,DX,[PERIOD],[WAVELENGTH]) <br>% Split S(SPACE,TIME) into: DW/ST/UP-waves<br>%<br>% Split the signal S(SPACE,TIME) into its downward/upward space <br>% propagating and stationnary components via a 2D Fourier decomposition.<br>% S is a (SPACE,TIME) matrix. We eventually proceed to a space and/or <br>% time filtering. DT,DX are temporal and space step. PERIOD=[Tmin Tmax] <br>% and WAVELENGTH=[Xmin Xmax] are time and space band-pass filters <br>% specifications.<br>%<br>% Sf(DIRspec,SPACE,TIME) are each components with:<br>%    DIRspec=1 -> Downward part (decreasing space axis)<br>%    DIRspec=2 -> Stationnary part<br>%    DIRspec=3 -> Upward part (increasing space axis)<br>%<br>% A(DIRspec,SPACE/2,TIME/2) are propagative Power Spectral Density with:<br>%    DIRspec=1 -> Downward part (decreasing space axis)<br>%    DIRspec=2 -> Upward part (increasing space axis)<br>% Rq: SPACE and TIME must be even<br>%<br>% FX and FT are space and time wavenumbers:<br>%    FX -> Space wavenumber, ie: 2*pi*[0:SPACE/2-1]/DX/SPACE<br>%    FT -> Time wavenumber, ie: 2*pi*[0:TIME/2-1]/DT/TIME<br>% <br>% PERIOD=[Tmin Tmax] is in real time axis unit<br>% WAVELENGTH=[Xmin Xmax] is in real space axis unit<br>%<br>% Rq: <br>%   - PERIOD/WAVELENGTH=[0 Inf] produce no filtering<br>%   - Tmin or Xmin = 0   will produce high-pass filtering<br>%   - Tmax or Xmax = Inf will produce low-pass filtering<br>%<br>% Ref: Park et al (2004) GRL V31<br>%      Park (1990) C. R. Acad. Sci. Paris 310(II) p919-926<br>%      J. Le Sommer, personnal communication<br>%================================================================<br>%<br>%  Guillaume MAZE<br>%  Last reviewed:<br>%                 Nathalie Daniault - September 2004<br>%                 -> Vectorisation of recomposition loop<br>%                 Nathalie Daniault - August 2004<br>%                 -> Change phases definition and <br>%                    adapt to a signal with non nul mean<br>%
             </td>
          </tr>
          <tr>
             <td style="border: 1px solid rgb(170, 170, 170); padding: 5px; background-color: #ddd; font-size:12px;">
                  Last update: 2011 March 04, 17:46<br>
                  Created by Guillaume Maze<br>
                  More informations at: <a href="http://codes.guillaumemaze.org/matlab">codes.guillaumemaze.org/matlab</a><br>
             </td>
          </tr>
     </tbody></table></center>
</body>
</html>
